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Classical chaos is often characterized as exponential divergence of nearby trajectories. In many 
interesting cases these trajectories can be identified with geodesic curves. We define here the entropy 
by 5 = Inx(x) with x(®) being the distance between two nearby geodesics. We derive an equation 
for the entropy which by transformation to a Ricatti-type equation becomes similar to the Jacobi 
equation. We further show that the geodesic equation for a null geodesic in a double warped space 
time leads to the same entropy equation. By applying a Robertson-Walker metric for a flat three- 
dimensional Euclidian space expanding as a function of time, we again reach the entropy equation 
stressing the connection between the chosen entropy measure and time. We finally turn to the 
Raychaudhuri equation for expansion, which also is a Ricatti equation similar to the transformed 
entropy equation. Those Ricatti- type equations have solutions of the same form as the Jacobi 
equation. The Raychaudhuri equation can be transformed to a harmonic oscillator equation, and 
it has been shown that the geodesic deviation equation of Jacobi is essentially equivalent to that 
of a harmonic oscillator. The Raychaudhuri equations are strong geometrical tools in the study of 
General Relativity and Cosmology. We suggest a refined entropy measure applicable in Cosmology 
and defined by the average deviation of the geodesics in a congruence. 


Classical chaos is generally defined as exponential di¬ 
vergence of nearby trajectories causing instability of the 
orbits with respect to initial conditions or quite simply 
as high sensitivity to initial conditions. The extent of di¬ 
vergence is quantified in terms of Lyapunov exponents 
measuring the mean rate of exponential separation of 
neighboring trajectories. 

The norm (for the Euclidean case, for example) 
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is a measure of the divergence of two neighboring trajec¬ 
tories, where Sxi is the z-th component of the displace¬ 
ment between two nearby trajectories at time t. The 
mean rate of exponential divergence is given by 
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The Kolmogorov entropy is related to the Lyapunov 
exponents. It gives a measure of the amount of informa- 
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tion lost or gained by the system as it evolves [9]. It can 
be computed from the Lyapunov exponent by 

hK= dfi, 

•>P ^.>0 

which is the sum of all positive Lyapunov exponents aver¬ 
aged over some region of the phase space P with measure 
dpi. 

One would naturally be interested in defining a mea¬ 
sure for stochasticity in regions with divergence. The 
function d(r) initially has an irregular behavior, and 
evolves into a form in which the limit as r —>■ oo of 

/ ( d{T) \ 
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converges to a value that depends on the initial condi¬ 
tions. Casartelli et al [3] argued that this quantity is 
deeply related to the Kolmogorov entropy and also ex¬ 
hibits strong sensitivity to the initial conditions. 

Benettin et al defined a similar entropy ^ and cal¬ 
culated a Kolmogorov-like entropy for the Henon-Heiles 
system. However we shall take a different route in this 
study. 

There are many interesting physically relevant exam¬ 
ples for which the trajectories can be put into correspon¬ 
dence with geodesic curves, for example, in problems in 
general relativity, and in the conformal map of Hamil¬ 
tonian potential models [6], with geodesic deviation de- 
scibed in terms of a Jacobi equation related to the cur¬ 
vature. In the following, we provide a relation between 
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the Jacobi equation, the entropy (as defined above) and 
the geodesic equation itself. 

Let there be given two nearby geodesics, m and n, and 
let T be the affine parameter on the geodesics. For a point 
X with parameter r on the geodesic m, one may define 
the geodesic deviation as the length of the shortest path 
from m to n. Let’s denote this geodesic deviation by 
x(r). The Jacobi equation states that [S] 

vl t1 

= -if(x(r))x(T) (1) 


where (j) = ((>(a;), a = a(t), and g is independent of t. 
Consider further a variation of the geodesic with Ja;* = 0 

Ss = J + aa'gijX^x^Stj dr = 0 

where ' = d/dt. Integrating by parts one gets + 

aa'gij0b‘x^)5tdT = 0 which implies 

+ aa'gijX^x^ = 0 , ( 10 ) 

For a null geodesic one has 


where K is the Gaussian curvature, and this simple form 
for the curvature is restricted to two-dimensional sys¬ 
tems. 

We employ the entropy defined by [1] 


5 = lnx('r) (2) 

From Equations Q-® we derive 

+ (<5 (t)) -t Js:(x(t)) = 0, (3) 

where ' = d/dr. One may transform Equation ([^ to a 
Ricatti-type equation by letting S{t) = X{t) 

X{T)+X{Tf + K{x{T))=0. (4) 


<ij.z z • 2 • 7 

t = a gijX x-’. 


Substituting into the geodesic equation (101 above leads 
to 


t + —r = 0. 

a 


( 11 ) 


We have achieved the equation which formally is the same 
as the entropy equation (^. 

Taking into account that the universe is evolving in 
time we study the entropy S{t) = Inx(T) in a four¬ 
dimensional cosmological spacetime with a time depen¬ 
dent metric. It is in fact a special case of a Robertson- 
Walker metric for a universe for which the space for a 
fixed time is a fiat three-dimensional Euclidian space ex¬ 
panding as a function of time The metric of the 

model is given by 


The general Ricatti equation has the form 


= -df + a{tf{dx^ + dy'^ + dz^). (12) 


X(T)=qo (t) -h qi{ t)X + q 2 {T)X'^. (5) 

The solution to equation ([^isX = —M/g 2 'a with u being 
the solution to the equation 


il — Til + Ru = 0 


( 6 ) 


The Christoffel symbols for the time components fi = 0 
are given by [5] 

r!^o = r°o = To* = 0, rl = a{t)dm,. 

By inserting these into the geodesic equation § we ob¬ 
tain 


In Equation § R = 9270 and T = 9i -f 92 / 92 ■ We there¬ 
fore obtain the equation 


d^a;° 


dr^ 


a(t)d{t)Sij 


dec* dx^ 
dr dr 


(13) 


u + Ku = 0, (7) 

This equation is the Jacobi equation in two dimensions. 
For a fiat space K — 0 and equation ([^ takes he form 


The Christoffel symbols for the spatial components (/i 
0) are 


pi _ pi _ ri 
^ jk ^00 


_ -pt _ 

^ Oj - 



S{r) + {S{r)r=0, ( 8 ) 

This equation has some resemblance with the geodesic 
equation 


dr^ dr dr 


(9) 


in particular if F))^ vanishes for p not equal to a. 

Consider, in particular, the geodesic equation for a null 
geodesic in a double-warped spacetime 


ds^ = —(jy^dt^ -\- a^gijdx^dx^ 


and the spatial part of the geodesic equation takes the 
form 


d^x* d(t) dx* dx^ 

-1—- 

dr^ a{t) ^ dr dr 


(14) 


^ The model used here provides a simple illustration of the similar¬ 
ity between the geodesic equation and the entropy equation 
which is the main intent of this study. Reference also splits 
the geodesic equation into a time part and a space part but uses 
a different technique with the purpose of obtaining the cosmo¬ 
logical redshift. 
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Equations (131 and (141 constitute the splitting of 
the geodesic equation into the timelike and spacelike 
parts [5]. 

For particles moving freely under purely gravitational 
forces one can find a freely falling coordinate system with 
the motion being a straight line in space time 


cPx^ 

dr'^ 


= 0 . 


(15) 


where r is the affine parameter along the geodesic, and 
Rab is the Ricci tensor of the metric [7], 0 = is 

the expansion, ~ the shear, and Wy = 

the twist. Round brackets represent symmetriza- 
tion, and square brackets represent anti-symmetrization. 

For completeness, and because it is very simple, we 
carry out the derivation of this equation explicitly. De¬ 
noting the covariant derivative by the geodesic 
equation is: 


Here t is the proper time 


= 0 , 


dr^ = iriapdx°' dx^ 


(16) and because = constant it follows that also: 


For massless particles the RHS of equation (16) van¬ 
ishes [2], and we may use u = cc® as the parameter 
instead of r. Photons follow null-geodesics and we re¬ 
strict ourselves to paths along the x-axis, i.e. x^{a) = 
{<(ct), a;((T), 0,0}. With the metric given by (12) and 
ds^ = 0 we obtain 


— dt^ + a{t)‘^dx'^ = 0. (17) 


This leads to the equation 


dx 1 dt 
da a{t) da 


(18) 


By solving for dt/da and inserting the null-condition (17) 
into the time component for the geodesic equation [4] we 
finally achieve the equation 


d'^t a{t) f 

da^ a{t) \dcr/ 


(19) 


This equation is formally identical to the entropy equa¬ 
tion 

It is noteworthy that a resemblance between the 
geodesic equation and the entropy equation is obtained 
by inserting the null condition into the time part of the 
geodesic equation and not the spatial part, which un¬ 
derlines the connection between the present definition of 
entropy with time rather than space. 

The definition of entropy as defined by equation ([^ 
has its origin in the geodesic deviation equation describ¬ 
ing the behavior of a one-parameter family of nearby 
geodesics and is, as remarked, in the present form re¬ 
stricted to systems of at most two dimensions. For higher 
dimensional systems one needs more refined tools to de¬ 
scribe the behavior of a bundle of geodesics, the so-called 
congruence. We now argue that the Raychaudhuri equa¬ 
tion may provide such tools in dimension 4. In a forth¬ 
coming study we shall show examples of entropy defined 
by the average deviation of the geodesics in a congruence. 

Let be the tangent vector field to a geodesic flow, 
and be the metric on the subspace perpendicular to 

The Raychaudhuri equation is 


de 

dr 



aija^^ + UJijU}’’^ 


Rij^ ^j, 


( 20 ) 


= 0. 

Without loss of generality we assume that = — 1. 
The metric on the spacelike subspace perpendicular to ^ 
is then 


dij — T ' 

We now decompose the derivative ^ij of ^ into three com¬ 
ponents: 


0 = 

^ + O'*) - ^0h^j, 

~ 2 ~ 0 *) ■ 


We note that the expansion 9 measures the logarithmic 
derivative of the volume element in the space perpen¬ 
dicular to the shear aij measures the non-conformal 
part of the defomation of the metric h, and the twist 
measures the entangling of the geodesic trajectories, i.e. 
the obstruction to ^ being hypersurface-orthogonal. The 
expansion 0 in (20) corresponds to .5 in (|^ and can be 
taken as the derivative of the entropy. Equation ([^ is 
the 2-dimensional version of (20). 

We can decompose: 


^ij — 2 0hij -)- aij iOij , 


and note that these three components are mutually or¬ 
thogonal: 

^ehija'-^ = = CTijO;*'* = 0. 

The first expression on the left vanishes because a is 
traceless, the second and the third vanish because hij 
and aij are symmetric, while uJij is anti-symmetric. 
Taking a derivative of 0 along we find 

0 = V 5 


where for simplicity we have denoted Vk^ij = ^ijk^ From 
the definition of the Riemannian tensor, we have 


^ijk ^ikj — 


p tm 

^jkims 5 
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hence we get 


0=- Rkm^r, 

where Rkm = g^^Rjkim are the component of the Ricci 
tensor. Also 




ij ck , 


Substituting the decomposition of ^ and using the or¬ 
thogonality relations, we obtain 

hik—ck^^ ckik+uJiktOik = —-^9^—a^+uP'. 


Substituting back into the equation for 9 we obtain (201. 


Consider now, the Einstein equations R^^, — \Rg^y = 
SttGE^^. Taking the trace we get R = —8nGT, hence 
substituting back into the Einstein equations we obtain 
Rk-iy = 8T^G{Tfj,i, - \Tg^^) and therefore R^^U^U'' = 
8'kG{T^i, — \Tg^i,)U^U''. Most known physical mat¬ 
ter fields satisfy the Strong Energy Condition (SEC), 
which states that for all time-like vectors C/, the inequal¬ 
ity > \Tg^„U^U'' holds. It follows, when the 

SEC holds, the term is always positive. Fur¬ 

thermore, note that the shear and the rotation are spatial 
vectors and consequently > 0, and > 0. 

As mentioned above, is zero if and only if the con¬ 
gruence is hypersurface-orthogonal. If that is satisfied 
the Raychaudhuri equation simplifies to the form 


+ a^ = -R^^u^ur 


( 21 ) 


In order for the LHS to be negative it must fulfill the 
condition dBjdr < —which finally leads to the in¬ 
equality: 


1 11 
- >- 1 - T 

9{t) - 00 3 


( 22 ) 


One concludes that any initially converging 
hypersurface-orthogonal congruence must continue 
to converge and within the finite proper time r < — 
will lead to crossing of geodesics (a caustic) which means 
that matter obeying the SEC cannot cause geodesic 
deviation but will increase the rate of convergence in 
accordance with the fact that the SEC causes gravitation 
to be attractive [5]. The aim to define the entropy by 
the average convergence/divergence of the geodesics in 
a congruence will be tantamount to establish that the 
SEC will cause initially decreasing entropy to continue 
to decrease. 

The Raychaudhuri equation for the expansion is a first- 
order nonlinear Ricatti equation and hence of the same 
type as equation (4) for which the solution. Equation 0. 
has the same form as the Jacobi equation. 


If we set 9 = 3F'/F the Raychaudhuri equation is 
transformed to 

^2 p 1 

— + - + a^-cu^)F = 0, (23) 

which is a harmonic oscillator equation. As pointed out 
above, 9 may be identified with the derivative of the en¬ 
tropy, so that according to (2) for the entropy S — In F 
here, F may be identified with an effective geodesic de¬ 
viation. 

We recently proved [H] that the geodesic deviation 
equation of Jacobi is essentially equivalent to that of 
a harmonic oscillator. The expansion 9 is the rate of 
growth of the cross-sectional area orthogonal to the bun¬ 
dle of geodesics. Increase/decrease of this area is the 
same as the divergence/convergence of the geodesics. 
The average growth of the cross-sectional area is com¬ 
patible with the average geodesic deviation. 

Kar and Sengupta have shown that the condition 
for geodesic convergence is the existence of zeroes in F at 
finite values of the affine parameter, and they argue that 
convergence occurs if 

+ a^-uj'^>0. (24) 


i.e. the shear accelerates convergence and the rotation 
obstructs convergence. 

Since shear transforms circles to ellipses, we compared 
the mean distance between uniformly distributed pairs 
of independent points inside a circle to that inside an 
ellipse of the same area, and found that it is smaller in 
the circle. The mean d of the distance d between pairs of 
points in a planar region fl of area tt was computed by: 

d=\ [ [ \x - y\dAj;dAy. (25) 

^ Jn Jn 


The result is graphed against the eccentricity e in Fig¬ 
ure [l] For comparison, we also computed the same quan¬ 
tity for rectangles of ’eccentricity’ e and area tt, where 
by similarity with the definition for an ellipse, we de¬ 
fined t he eccent ricity of a rectangle with sides a > 6 as 
e = conjecture that the mean distance 

between pairs of points inside any plane domain of area 
TT is smallest for a circle, i.e. the circle is the unique min- 
imizer of (25) among all planar regions of area tt. 

To the best of our knowledge, this conjecture has not 
been addressed before. Nevertheless, it might have im¬ 
portant implications. The evolution from an infinitesimal 
circular cross-section orthogonal to the flow lines to an 
elliptical one of same area is brought about by shear. 
Moving the cross-section along the flow does not change 
the number of geodesics. However, due to increase in the 
mean-distance between the geodesics when transforming 
from a circular to an elliptical cross-sections, there is a di¬ 
verging tendency of the geodesics moving along the flow. 
That implies, according to our proposed definition of en¬ 
tropy as the mean distance between geodesics in a bun¬ 
dle, that the evolution in the presence of shear exhibits 
an increase of entropy. 
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